Doping dependence of the electron-doped cuprate superconductors from the 
antiferromagnetic properties of the Hubbard model 
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Within the Kothar-Ruckenstein slave-boson approach, we have studied the antiferromagnetic (AF) 
properties for the t-t'-t"-U model applied to electron-doped cuprate superconductors. It is found 
that, due to the inclusion of quantum fluctuations, the AF order decays with increasing doping 
much faster than obtained in the Hartree-Fock theory. Under an intermediate constant U the 
calculated doping evolution of the spectral intensity in the AF state has satisfactorily reproduced 
the experimental results, without need of a strongly doping-dependent U as argued earlier. This 
may reconcile a discrepancy in recent studies on photoemission and optical conductivity. 
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I. INTRODUCTION 

The intriguing Fermi surface (FS) evolution with 
doping in electron-doped cuprate Nd2_xCea;Cu04 re- 
vealed by angle-resolved photoemission spectroscopy 
(ARPES) measurementsi*^ has attracted much atten- 
tion recentlyi^i'^i^i^i'^i^i^i^°i"i^^ It was observed that at low 
doping a small FS pocket appears around (tt, 0). Upon 
increased doping a new pocket begins to emerge around 
{tt/2, it/2) and eventually at optimal doping x = 0.15 the 
several FS pieces connect to form a large curve around 
(tt, tt). By use of the t-t'-t"-U model and Hartree-Fock 
(HF) mean-field treatment, Kusko et aim' have studied 
the FS in consideration of the antiferromagnetic (AF) 
order and reached remarkable agreement with experi- 
ments. However, in their work the on-site U is treated 
as a doping-dependent effective parameter, specifically 
U — 6t at doping x = and 3.1t at a; = 0.15. Simi- 
larly a doping-dependent U was argued by others based 
on numerical calculationsi^ Alternatively, the strong- 
coupling t-t'-t"-J model was adopted by some of us to 
construct the FS.^ Without tuning parameters we have 
also obtained consistent results with ARPES data. 

Very recently, a systematic analysis of the optical con- 
ductivity for electron-doped cuprates was undertaken by 
Millis et al}- It suggests that (i) the electron-doped ma- 
terials are approximately as strongly correlated as the 
hole-doped ones and (ii) the U value is not strongly dop- 
ing dependent within the electron-doped family. Both 
of these implications are in apparent disagreement with 
the theoretical studies by Kusko et al.^ and in favor of 
the strong-coupling model used by us^ On the other 
hand, our result at optimal doping from the t-t'-t"-J 
modefi has not been reconciled by exact diagonalization 
calculations lii Thus it is still open whether the strong- 
coupling model, as implied by point (i) above, is available 
to explain the ARPES data at optimal doping. 

Here our attention is focused on the issue relevant to 
point (ii) above, i.e., for the t-t'-t"-U model whether the 
strong doping dependence of U as argued by Kusko et 
almi^ is really necessary for the interpretation of ARPES 



data. Theoretically it is not convincing to have a strongly 
doping-dependent U since it is a local interaction within 
an atom and difficult to be changed by adding carri- 
ers. To understand why a strongly doping-dependent U 
is needed by Kusko et at, we come to some details of 
their work. Within the HF mean-field theory, two AF 
energy bands are derived, with the minimum of the up- 
per band around (tt, 0) and the maximum of the lower 
one around (7r/2,7r/2). At low electron-doping, the up- 
per band is crossed by the Fermi level, leading to a small 
FS pocket around (7r,0). With increasing doping, the 
AF order weakens and the lower band shifts towards the 
Fermi level. The essential observation is that the AF or- 
der decays with doping too slowly under a constant U, 
and consequently the lower band is still too far away from 
the Fermi level even at optimal doping to contribute any 
spectral intensity around (7r/2,7r/2). Therefore an effec- 
tive U{x), which decreases with increasing x, is enforced 
to expedite the decline of the AF order so that the lower 
band approaches the Fermi level rapidly leading to the 
consistency with ARPES results. 



Actually, as we know, the HF theory overestimates 
the AF order due to the ignorance of fluctuations. So 
a tempting idea is that the real AF order should decay 
with doping in a (much) quicker way than obtained in the 
HF theory. Then it is interesting to ask how much the 
HF results will be improved if the fluctuations are taken 
into account and whether the ARPES data can already 
be explained with a doping-independent U. In this paper 
we try to answer these questions by studying the t-t'-t"- 
U model with the Kotliar-Ruckenstein slave-boson (SB) 
approach>i4 The method well considers the AF fluctua- 
tions for a wide range of interactionsJ^-- and has been 
used to study the stripes in hole-doped cuprates.^^ We 
have found that, compared to the HF results, the AF 
state derived from the SB approach has a much lower 
energy and the order parameter decreases much faster 
with doping. Consequently, the ARPES results can be 
reproduced under a constant U. 
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II. FORMALISM 

We start with the t-t'-t"-U model on a square lattice 
which reads 

{ij)a {ij)2<T 

't" J2 (cl^J- + ^^-c-) + ^ "^T«a : (1) 

where (), ()2, Os represent the nearest neighbor (n.n.), 
second n.n., and third n.n. sites, respectively, and the 
rest of the notation is standard. Throughout the work t is 
taken as the energy unit and typical parameters t' = —0.3 
and t" = 0.2 are adopted. 

In the spirit of the Kotliar-Ruckenstein SB approach, 
four auxihary bosons e[^\ p^^J {a ='\, |), d-^'' are intro- 
duced at each site to label the four different states, which 
can be empty, singly occupied by an electron with spin up 
or down, or doubly occupied. The advantage of introduc- 
ing bosons is that the interaction term can be linearized, 
i.e., the product of density operators n.i|n.i| is replaced 
by the occupation number operator for double occu- 
pancy. As the expense, the hopping term becomes com- 
plicated because any hopping process of electrons must 
be accompanied by the transitions of slave bosons. Ex- 
plicitly, if an electron (with spin a) hops from site i to 
j, the slave bosons must change simultaneously at both 



sites i and j. For example, at site i, associated with the 
annihilation of the electron the bosonic state will tran- 
sit either from the singly occupied one with spin a to 
the empty one or from the doubly occupied one to the 
singly occupied one with spin a. Namely, the transition 
of the bosons at site i can be described by the operator 
Zia = e\pia+p\gdi. To eliminate the unphysical states in 
the enlarged (fcrmionic and bosonic) Hilbert space, the 
following constraints have to be imposed at each site: 

eIe» + 5IpLP- + 4««» = 1 . (2) 

pLpit + d\d,, = c\^c,„ . (3) 

The first constraint states that a site is either empty, 
singly occupied, or doubly occupied. And the second one 
guarantees that, when an electron with spin a locates at 
site z, this site is either singly occupied with spin a or 
doubly occupied. 

In order to study the AF order, we divide the lat- 
tice into two sublattices A and B. For sublattice L 
{L = A,B), we introduce a set of bosons e[\^ , p^^j^^ {a =] 

, I), and Lagrange multiphers A.J^.' Kl which are 

associated with the constraints (O and ©, respectively. 
Thus the original Hamiltonian, with Lagrange multiplier 
terms added, can be rewritten, in terms of fermionic op- 
erators a and h corresponding to sublattice A and B re- 
spectively, as follows 



H 



-t (4Aa(4abj<yZjBa +ii-C.) ~t ^ (4i3<TCaj<T^jVl<T + h.c.) - i' ^ (z|^^a|^aj^ZjAa + h.c.) 

{ij){T,iGA {ij)a,i^B ('tj)2cr,iGA 

(ij>20',ie-B {ij)3a,i£A {ij)3a,ieB 



ieA 



ieB 



ieA 



ieA,a 



" 4b(^^B - Y^BaPiBa " 4gd,B) + Y ^iB^a^io ^ p\B<yPiBa ~ 4Bd^B) 



ieB,a 

r 



(4) 



where the operators 



elxPtLa + plLgd.L 



(1 - el^e^L ~ p\lsPiLs){1 - d\^d,L - p\LaP^La) 



have been renormalizedM to ensure the exact result at 
[/ = even in mean-field treatment as adopted sub- 
sequently. With mean-field approximation, the bosons 
are replaced by c numbers and assumed to be site- 
independent on each sublattice, i.e., (e-l') — Cl, {PiLa) ^ 
PLa, (d-i) = dL (and correspondingly = ZLa). 



At the same time, the constraints are softened to be 
satisfied only on the average on each sublattice, i.e., 
Aj^ \\, Af^ — > A£. This treatment is equivalent to 
making a saddle-point approximation in the path-integral 
formulation. Now the Hamiltonian becomes quadratic in 
terms of fermionic operators and can be easily diagonal- 
ized in momentum space. Finally we have 

H = Y^^kA.^k. + ElhlM + Eo , (5) 



where the energy bands read 



3 



(6) 



with 

£k = — 2t(cOS fca; + COS fcj,) , 

e'l^ = — 4t' COS fca; COS ky — 2t"(cos 2/?^: + cos 2ky) , 
and the constant 



Eo = (A^/2) 



C/(4 + 4) - Aji(ei + + 4 - 1) 
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1) 



2 

Bcr 



4) 



Above TV is the total number of lattice sites, k is re- 
stricted to the magnetic Brillouin zone (BZ), and the 
operators d and b are related to a and b by unitary trans- 
formations. 

The grand canonical thermodynamic potential is given 
by [/3 = l/iksT) with T: temperature] 



1 



(7) 



Here fj, is the chemical potential decided by — {d^ /9fJ.) = 
{Nf. is the total number of electrons). All the pa- 
rameters are determined to give the lowest free energy 
F = n + iVe/i. The calculation will be simplified if we 
search for the AF solution which satisfies the relations: 

SA — SB, dA ^ ds, PAa — PBS, — ^B^ ^A — ^B- 

(Correspondingly the energy bands E^^ become spin de- 
generate, i.e., i?^ — Ef.) The staggered magneti- 
zation, which characterizes the AF order, is given by 
m = {p% -4i)/2 = (4i -4t)/2- 



III. RESULTS 



In the following we are limited to T = only at 
which the AF long-range order is possible for the two- 
dimensional model. 

To look at the differences between SB and HF results, 
we first consider the half-filled case, i.e., n = N^/N = 1. 
As is well known, the pure t-U (with t' = t" = 0) model 
takes on the AF order for arbitrarily small U > due 
to the perfect nesting effect. After inclusion of the long- 
range hoppings, the AF order will be frustrated and a 
critical U — Uc is required to stabilize it. The relations 
between Uc and t' , t" from both SB and HF approaches 
are compared in Fig. ^ While Fig. ^a) shows Uc vs 
t' when t" = which has been given previously, we 
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FIG. 1: (Color online) The critical Uc for the AF order to be 
stabilized at half filling; (a) Uc vs —t' when t" — (b) Uc 
vs t" at fixed t'/t = -0.25 [blue (black) lines] and -0.3 [red 
(dark gray) lines]. All solid lines are from SB approach and 
dashed ones from HF theory. 



further present Uc vs t" at fixed t' in Fig.^b). It is clear 
that Uc increases with both t' and t". Moreover, the value 
of Uc at fixed t' and t" obtained from SB approach is 
always larger than that from HF theory. This is expected 
in view that the HF theory overestimates the AF order. 
Once the fiuctuations are taken into account, the critical 
value Uc will be largely raised as seen in the SB approach. 
We notice that for parameters t' — —0.25 and t" = 0.1 
as adopted by Kusko et al^ the accurate critical value 
is Uc = 4.23 from SB approach rather than 3.05 from 
HF theory. Thus those values for doping-dependent U 
selected by them^ i.e., U = 3.5(3.1) at x = 0.1(0.15) 
are actually unreasonable, which are already too small 
to stabilize the AF order. 

Then we come to see how the AF order declines with 
increasing doping x{— n — 1) under fixed U. For t' = 
—0.3 and t" = 0.2, the critical value at half filling is 
Uc = 5.64(SB)/3.81(HF). We choose f7 = 6.3 > Uc and 
solve the self-consistent equations. The staggered mag- 
netization m as a function of x is shown in Fig. Ola), 
where the red (dark gray) and black lines are obtained 
from SB and HF approaches, respectively. It is seen that 
at half filling the AF order derived from SB approach 
becomes weak relative to that from HF theory. And it is 
even weaker with increasing doping because m decreases 
much faster in SB approach than in HF theory. In addi- 
tion, we have found that the AF state becomes energet- 
ically unfavorable at a; ~ 0.157 in SB approach, leading 
to a first-order transition as shown by the dashed line 
in Fig. |2Ia). This first-order termination of the AF or- 
der, compared to the claim of a quantum critical point 
seems to be closer to the experimental observation where 
a steep disappearance of the AF phase was found. Also 
the doping range that the AF order survives is quantita- 
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FIG. 2: (Color online) Staggered magnetization (a) and 
ground state energy (b) vs electron doping for parameters 
t' = -0.3, t" = 0.2 and U = 6.3. The red (dark gray) 
lines are from SB approach and black ones from HF theory. 
The dashed line in panel (a) shows that m jumps to zero at 
X ~ 0.157 in SB approach, and correspondingly the dashed 
one in panel (b) gives the energy at m = 0. 

tively consistent with experiment)^ In contrast, m keeps 
finite up to a; ^ 0.4 and goes to zero continuously in HF 
theory (not shown) . The corresponding ground state en- 
ergies from both approaches are compared in Fig. |2Ib). 
The energy derived from SB approach is notably lower 
than that from HF theory and the energy difference in- 
creases with increasing doping. This indicates that the 
study of the AF property is largely improved within the 
SB approach. 

Subsequently the energy bands — fj, are plotted in 
Fig. 121 with red (dark gray) lines from SB approach and 
black ones from HF theory. (The definitions of in 
the magnetic BZ have been analytically extended to the 
original BZ.) While the upper bands derived from two 
approaches look similar at each doping, the lower bands 
take on different behaviors. Within HF theory, the lower 
band shifts very slowly towards the Fermi level with in- 
creasing doping and is still far away from it at x = 0.15. 
On the other hand, within SB approach, the lower band 
lifts towards the Fermi level as a whole. At a; = 0.1 it 
is already rather close to the latter and may contribute 
the spectral intensity within the experimental resolution. 
At optimal doping, it is nearly crossed by the Fermi level 
around (7r/2, n/2). 

In order to directly compare with the ARPES data, 
we have calculated the spectral function A{k,uj) within 
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FIG. 3: (Color online) The energy bands at different dopings. 
In each panel, the red (dark gray) and black lines are from SB 
and HF approaches, respectively. The Fermi energy is fixed 
at zero. 

SB approach. The density plots for integration oi A{k, uj) 
times the Fermi function over an energy interval [—40, 20] 
meV (same as that adopted in ARPES experimentsi) 
around the Fermi level are shown in Fig. ^ by the top 
row. At low doping x = 0.04 a small FS pocket forms 
around (tt, 0) [and equivalently (0, tt)]. At doping x = 0.1 
the spectral intensity begins to appear around (7r/2, 7r/2), 
and becomes strong at optimal doping. At the same time, 
the Fermi patch around (tt, 0) deforms with increasing 
doping, i.e., half of the FS around (7r,0) loses much of 
its intensity. The theoretical results agree well with the 
ARPES datai For comparison, the corresponding plots 
from HF theory are shown by the bottom row in Fig. 01 
which, however, do not exhibit the main experimental 
features. 



IV. DISCUSSIONS AND CONCLUSION 

In the above sections, after suitable consideration of 
the AF fluctuations we have reached consistent results 
with ARPES measurements under a doping-independcnt 
constant U, in contrast with the previous argument for a 
doping-dependent U j '^i'^i^i^ The agreement is satisfactory 
in view that all input parameters except U are extracted 
from experiments. We do not exclude that a tunable U 
with doping will help lead to a perfect comparison be- 
tween theory and experiment. Even in this case, how- 
ever, not a strong doping-dependcnce of U as given by 
Kusko et al. (~ 50% variation from x = to a; = 0.15)^ 
is needed. This may reconcile one of the discrepancies in 
recent explanations of both ARPES and optical conduc- 
tivity mentioned before. 

As for the magnitude of U, an intermediate value 6.3t 
is adopted here, which is somewhat too small at x = to 
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FIG. 4; (Color online) Density plots for integration of 
A{k, ij) times the Fermi function over an energy interval 
[-0.12, 0.06]i ~ [-40, 20] meV around the Fermi level, at dop- 
ing X = 0.04, 0.1 and 0.15 (from left to right). The top row is 
from the SB approach and the bottom row from HF theory, 
under a doping-independent constant U = 6.3 (with t' = —0.3 
and t" = 0.2). The Lorentzian broadening is 0.154. 



the AF fluctuations are fully taken into account, an even 
quicker decreasing of m with doping than displayed in 
Fig.EIa) will be the case. Then a larger U than adopted 
here, towards strong-coupling regime, could be chosen 
and fixed throughout the doping range, leading to the 
similar intensity maps shown in the top row of Fig. 01 
The evaluation of the biggest U available to explain the 
ARPES data as well as the possibility to cross over to 
the strong-coupling model need further work^ 

In conclusion, we have studied the AF properties for 
the electron-doped t-t'-t"-U model by use of the Kotliar- 
Ruckenstein SB approach. With inclusion of quantum 
fluctuations the AF order declines with increasing doping 
much faster than obtained in the HF theory. The calcu- 
lated doping evolution of the spectral intensity under a 
constant U is in good agreement with ARPES measure- 
ments on electron-doped cuprates. A strongly doping- 
dependent U, which is difficult to justify itself, is found 
to be actually unnecessary. 
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